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The mineral fuels used recently, i.e., oil and gas, will be soon exploited out. This paper presents an idea of 
the engine where any fuel or solar heat can be used as a source of energy. The proposed model is an 
externally heated, 2-stroke, valve engine (EHVE). This is a piston-type engine, entirely different from the 
well-known Stirling one, which is the best known example of such a solution. It works in a closed Joule 
cycle and is designed to produce a moderate amount of energy. The engine is composed of typical parts 
met in piston designs: an expander, a compressor, a heater, a cooler and, additionally, two recirculation 
blowers, which consume a small amount of produced power. An additional advantage is its working 
medium, which may be simply atmospheric air and the engine has a conventional crankshaft and an oil 
lubrication system. It has already been proven that operation of the EHVE is possible with satisfactory 
power and efficiency at the output. Comparisons of the EHVE action with and without recirculation 
blowers are performed. 

Crown Copyright © 2010 Published by Elsevier Ltd. All rights reserved. 


1. Introduction 

Our earlier publications [ 1 -4] have been devoted to an externally 
heated two-stroke piston air engine (EHVE), entirely different from 
the well-known Stirling engine, [9-11 ]. This paper is also devoted to 
a 2-stroke type of the EHVE and the major difference between it and 
the previous propositions refers to the magnitude of heat 
exchangers. They were small if compared to the volume of the 
working cylinder ( Vu/V E < 1) and now are much bigger - Vu/V E = 3-5. 
Additionally, two centrifugal blowers (not present in [1-4]) and 
a physical model of the complete engine have been added. The 
blower actions are included in the calculations and are assumed to 
consume rather a small amount of power, if compared to that 
produced by the engine. The power consumed by the centrifugal 
blowers increases monotonously with their dimensions and we 
assume that it is limited to 10% of the power produced by the engine, 
but the flow in the exchangers is turbulent. A kind of optimization of 
power consumed by blowers can be found in Fig. 9. 

An important advantage of the proposition is use of the easiest 
available working medium - atmospheric air and a possibility to 
apply any source of heat, i.e., any fuel including solar heat or nuclear 
power. As most mineral fuels used recently, i.e., oil and gas, will be 
exploited out soon, the idea seems to be worth further develop¬ 
ment. The EHVE is designated to produce a moderate amount of 
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energy on the basis of an external source of heat. It consists of a few 
essential parts: an expander, a compressor, a heater, a cooler, two 
recirculation blowers and works in a closed Joule cycle. The novelty 
proposed in the paper is an application of two blowers working in 
loops of high and low pressures, which are linked to the engine 
heater and the cooler. The air circulation has to be ensured in the 
exchangers when the valves of the EHVE are closed and the blowers 
exist. The engine [8], its expander and compressor, are similar in 
design to those presented in [1-4], though our solutions have been 
developed independently. Behaviour of heat exchangers used in the 
experiment [8] was not described, though their action is crucial for 
the engine operation. 

Brzeski and Kazimierski [5] present results of the experimental 
investigation of the EHVE prototype realizing the idea described in 
[1-4]. Two heaters of small volumes were employed. The heater 
contacts the expander and the compressor only during 1 \2 period 
of the engine cycle. During the remaining time, isochoric heating of 
air enclosed in the heaters should take place. The experiments 
show, however, that such a type of heating is insufficient for proper 
operation of the EHVE. The working gas must be changed. Any gas 
which takes heat by radiation should be applied, and the wall 
temperature of the heater has to be high enough. Therefore, these 
changes in the present paper have been introduced. 

The issue of the 4-stroke design of the EHVE should be 
considered thoroughly. The problem is connected with the time 
when a valve is opening. In 2-stroke engines, the time corresponds 
to an angle of about 90°, but in the case of 4-stroke engines the 
camshaft rotates twice as slow. All valves in the 4-stroke engine 
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Fig. 1. Scheme of the transversal cross-section of the engine, E - expander, 
C - compressor, H - heater, CL - cooler, Bl, B2 - blowers, 1, 2, 3, 4 - valves. 

should be composed of two typical valves like those presented in 
[6] and such a solution is very difficult in practical realization. The 
experimental investigation of the 4-stroke EHVE was performed in 
2003 and the results were not encouraging enough and the main 
reason was a very complicated structure of the valves. The pre¬ 
sented return to the idea of the 2-stroke engine with a usual value 
of volumes of heat exchangers and blowers is a way of development 
of the EHVE when no radiation is taken into account. 

The paper is organised as follows. Section 2 presents a principle 
of the engine operation, a complete theoretical model is described 




Fig. 3. Scheme of air flow in the heater loop. 


in Section 3, results of computer simulation are presented in 
Section 4, and finally conclusions are shown in Section 5. 

2. Principle of the engine operation 

The described EHVE represents a 2-stroke-type engine. It is 
equipped with two cylinders, each of them playing a different role: 
an expander (E) and a compressor (C) situated between two heat 
exchangers of moderate volumes: a heater (H) and a cooler (Cl). The 
phase angle between compressor and expander pistons is equal to 
180°. Both the heat exchangers are additionally supplied with their 
internal blowers raising the rate of heat exchange. The engine 
works in a closed cycle of the Joule type and its working medium is 
simply air. The schemes of the EHVE are shown in Figs. 1 and 2. A 
scheme of flows of the working medium, when the blowers exist, is 
shown in Figs. 3 (for the heater) and Fig. 4 (for the cooler). Figs. 3 
and 4 are added for better understanding of the collector 
volumes. All valves of the engine are operated by a traditional 
camshaft. Additionally, the EHVE has a conventional crankshaft and 
an oil lubrication system. 

When the piston of the cylinder E starts to move down from its 
upper position, valve 1 opens and the working medium at high 
pressure and temperature values flows from the heater H to the 



Fig. 2. Scheme of the longitudinal cross-section of the engine. 


Fig. 4. Scheme of air flow in the cooler loop. 
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Fig. 5. Pressures in essential parts of the engine at 3000 rpm as a function of the 
crankshaft angle. 

expander E - Figs. 1 and 2. Valve 1 is doubled (it means, its action is 
realized with 2 identical valves) for better filling of the cylinder E. 
Both valves 1 close at the angle equal to oq c and further piston 
movement causes an expansion of the hot air in the locked volume 
of the cylinder E until the piston reaches its lowest position. At the 
same time the air locked in the volume of the cylinder C (all the 
valves closed, see Figs. 1 and 2) is compressed and prepared for 
flowing to the heater. 

Valve 3 opens at about 100° and the compressed air flows into 
the heater under high pressure - Figs. 1 and 2. In this design, two of 
valves 3, shown in Fig. 1, are employed for better transportation of 
the compressed air to the heater H. When the piston C reaches its 
upper position, the supply of the compressed air to the heater 
finishes (Fig. 1), and valves 3 close. Then, the pistons in the cylin¬ 
ders E and C start to move again. The piston in the cylinder E moves 
up and valve 2 opens, allowing the cooler to take the overworked 
air (Fig. 2). At the same time the piston in the cylinder C moves 
down, valve 4 opens and the cold air from the cooler flows to the 
cylinder C. Finally, the pistons reach their primary positions, valves 
2 and 4 close and the two-stroke engine cycle ends, realizing a full 
Joule cycle. 

The direction of opening of valves 1 and 3 (Figs. 1 and 2) is 
inverse than valves 2 and 4, which act in accordance with valves 
usually met in internal combustion engines. 



Fig. 6. Temperatures at 3000 rpm - explained at the end of Section 4 as a function of 
the crankshaft angle. 



Fig. 7. Mass flow rates of the engine valves for 3000 rpm as a function of the crank¬ 
shaft angle. 

Publications [1-4] and [8] present a compressor and an 
expander combined into a single device, which yields a compact 
cylinder set, however it is difficult in practical realization although 
Vc/V E = 1. Such a solution is also possible in the presented case, but 
it could be even more difficult in practical realization because in 
that case Vq/Ve < 1. Besides, in [3] a thermodynamic cycle of the 
previous version of the EHVE is described. Two self-acting valves (3 
and 4) were used in the above-mentioned propositions. 

3. Theoretical model 

The model presented below is based on the statement formu¬ 
lated in the introduction. All additional assumptions are quoted 
later in the paper in places they are introduced. We assume 
approximately that the engine cylinder (cooled) operates in such 
a way that a heat exchange inside it does not exist. 

The most important for cylinders are the time-dependent rela¬ 
tions and this has been taken into account, whereas a dependence 
on spatial variables has been neglected, so the model is described 
and simulated in the time domain only. 

The heat exchangers are situated between two collectors, see 
Figs. 3 and 4. The inlet and outlet collectors are excluded from the 
heat exchangers and it is assumed that no heat streams inside the 
collectors exist. The collectors play, at the same time, a role of 
settling chambers, which are important for the behaviour of 



volume [m 3 ] 


Fig. 8. p-V diagram for 3000 rpm. 
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installation. The accuracy of the approach is checked many times 
during the computer simulation in Section 4. 

Such a treatment of the parameters is of course a kind of 
approximation, but it allows us to use a moderate number of 
equations in case of a time-dependent model. These are non-linear 
ordinary, differential equations of mass and energy conservation, 
supplemented by some non-linear algebraic ones. The details of the 
theoretical model will be explained later in the paper. 

Firstly, let us consider geometry of cylinders and valves, then 
mass flow rates through the EHVE valves will be presented and, 
finally, essential differential equations for the main elements of the 
engine will be discussed. 


3.1. Geometry of cylinder volumes and the valve cross-section area 


The volume of the expander is described by: 


S E 


ft()E + 2~(1 ~ cos 


where: A E = 7rd?/4, and its derivative is: 

dVr SnCxJ . 

— - 7 — = A e -F— sin (ot. 

d t 2 

The volume of the compressor equals to: 


Vr = A ( 


Sr* 

hoc + -£(1 +cos tot) 


( 1 ) 


( 2 ) 


( 3 ) 


where: Ac = 7rdc/4, and its derivative is: 

d Vr Sr(x) . 

—= -A c -±—sin cot. 

dt 2 

The cross-section areas of the valves - Aj {j = 1, 2, 3, 4) 
described by the following formulae (where a = ot): 


( 4 ) 

are 


1) for aj t 0 < oi < aj t 0 + A af 


Fig. 9. Power (a) and efficiency (b) of the engine as a function of its rotational speed. 


blowers. The preliminary calculations show that pressures and 
temperatures in the collectors are almost constant (Figs. 5 and 6). 
Such a situation takes place when the rotation of the engine is 
constant and a volume of the collector is greater than a volume of 
the cylinder. This allows for considering the blowers as almost 
stationary devices. This is explained in more detail in Section 4 of 
the paper. The power consumed by the blowers will be determined. 

In both heat exchangers the spatial changes may be important and 
they were earlier calculated in [4], as space and time variables, 
because the temperature of the working air inside the heater in fact 
grows from the inlet to outlet collectors - Figs. 1 and 3. A simpler 
approach to the problem is also possible and it is applied in this paper. 

The heat exchangers can be described as follows. The proper 
heater and cooler are counter-current heat exchangers, see Fig. 1, 
where a temperature drop remains almost constant along their 
length, [12]. The heat exchangers are treated as time-dependent 
only, assuming that all important changes are the same along the 
overall length of the exchangers and considering each of them in 
a single plane only. The plane is situated at the outlet of the 
exchanger tubes. The approximate approach used here allows for 
applying maximal and minimal temperatures of the exchanger walls 
in the discussed case. The heat volumes delivered to and subtracted 
from the EHVE in the heat exchangers are estimated accordingly. 

The applied approach to the exchangers is sufficient for calcu¬ 
lations of the benefits resulting from using the blowers in the 


1 / CL — 0L\ n 

Aj( a ) = 2 max ( 1 - cos ^ 



2) for aj t 0 + A aj < a < a JjC — A af. 


Aj( 0 i) — -4j,max» (6) 

3) for aj' C - Aaj < a < aj tC : 


1 ( 0i - OLi c \ 

Aj( a ) = 2^fmax ( 1 - cos ^ ’ J 1 

and: 



^j.max — tt^(/ • hj • sin 45 . (8) 

All the data: aj t 0 , oij iC , Aaj , dj, hj are assumed in the computer 
program. The last term in Eq. (8) results from a typical design of the 
shape of valves. 


3.2. Mass flow rates through the engine valves 

Mass flow rates are calculated using Bernoulli and continuity 
equations for gas, where the losses are determined by coefficients 
Zj. The general form of this equation (for j = 1, 2, 3, 4) is: 
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riij = 
and, if: 


2x 


K — 1 


1 



( 9 ) 



( 10 ) 


• _ 2 
J v/RTb U+l 


/c+1 

2(k-1) 


Values of the coefficients equal to 0.85 are usually used in 
similar designs. These values are found in literature; see for 
example [14]. Pressures p a and pb and temperatures T a , 7b in Eqs. 
(9)—(11) are calculated as follows: 


for j = 1 

•f Phi > Pe = Pb 

if Pe > Phi : Pb 
for j = 2 

if Pe > Pci 2 : Pb 

if Pci 2 > Pe = Pb 
for j = 3 

if PC > PH3 : Pb 

if Ph3 >Pc-Pb 
for j = 4 

if PC > Pci4 : Pb 

if Pci4 > Pc ; Pb 


Phi, 

T b - 

= I'm, . 

Pa 

= Pe, 

Pe, 

Tb = 

Te, Pa 

— 

Phi, 

Pe, 

Tb = 

: Te- Pa 


PC 12 , 

PC 12 

, T b 

= 7 ci 2 > 

Pa 

= Pe 

Pc, 

T b = 

Tc, Pa 

= 

PH 3 , 

PH 3 , 

T b - 

= 7 h 3 > 

Pa 

= Pc, 

Pc, 

T b = 

: To Pa 

— 

PC 14 , 

Pci 4 

, T b 

= 7 ci 4 > 

Pa 

- Pc 


( 12 ) 


The mass streams rh H i and rh H 3 between parts of the heater will 
be discussed afterwards. 


If the equation of state for real gas (e.g., Peng-Robinson) is taken 
into account, the error for the extreme parameters of air (100 bar and 
1000 °C) is less than 2%. So, we remain with the quoted equation. 

In case of the compressor, we have similar dependencies. Its 
equation of mass conservation takes a form: 


1 ( . . dV c 

±m 3 ±m 4 - p c - 


dp c 

d t ~~ dt 

The signs of the mass streams m 3 and rh 4 are given below: 


(17) 


// 


( 11 ) ™3 


if Pc>PH 3 , sign 
if Ph 3 > pc, sign " +" 


m . if Pc > Pci4, sign " 
if Pci 4 >Pc, sign 


// 


// _| _// 


(18) 


The equation of energy conservation for the compressor i 
follows: 


is as 


dl c 

dt 


-T c (±m 3 ±m 4 ) +k 
Pq dV c 


f +^3^H3 
1 -tn 3 T c 


/ +rh 4 T CL 
+ l ~m 4 T c 



c v dt 


(19) 


Eqs. (15) and (19) assume that no heat exchange takes place in 
the cylinders of the EHVE. The signs of the mass stream terms m 3 
and m 4 are determined in Eq. (18). The pressure pc inside the 
compressor is calculated from the equation of state: 


Pc - PjcR- 


( 20 ) 


Let us consider the heater model now. It is divided into three 
parts of constant volumes, namely: 


3.3. Essential equations for EHVE parts 

Let us consider the expander first. In its volume V E , the equation 
of mass conservation has a form: 


dt - V E 


1 / , . dV E 

±m 1 ±m 2 - p E -^ 


(13) 


The signs of the mass streams rrq and rh 2 are determined as 
follows: 




m 2 


if 

Phi >Pe, 

sign 

if 

Pe >Phi, 

sign 

if 

Pe > Pci 2 , 

sign 

if 

Pci 2 > Pe, 

sign 


n _// 

//_// 
// i // 


(14) 


The equation of energy conservation for the expander takes 
a form: 


dip 


1 


dt p E V E _ 


T E (±m 1± m 2 ) +K ({+^ 


f + rr*2 Tci2 1 A Pe d^ E 
l -m 2 T E J 


'V 


dt 


(15) 


The signs of the mass stream terms rrq and rh 2 can be deter¬ 
mined from Eq. (14). The temperature in the parenthesis is chosen 
accordingly to the relations expressing the mass flows. Such 
a notation is repeated in the following equations of energy 
conservation. The pressure p E inside the expander is calculated 
from a typical equation of state for ideal gas: 

(16) 


1. inlet collector - V H 3 , 

2. proper heater (a counter-current type) - Vh, 

3. outlet collector - Vhi- 

The volume of the proper heater and its surface are determined 


by: 

V H = 

A h = it d H / H n H 


( 21 ) 


The working air flows from the compressor C to the inlet collector 
Vh 3 through valve 3 (Fig. 3). Then, it flows to the proper heater volume 
Vh, transfers to the outlet collector Vhi, and finally reaches the 
expander E volume through valve 1 - Fig. 3. The blower B1 transports 
the air from the outlet collector of the heater to its inlet. We assume 
that the temperature of the outlet collector is also transferred to the 
inlet one and the pressure drop between them is maintained by the 
blower Bl. The heater collectors are considered as settling chambers 
for the flow. The mass flow rates rh H i and m H3 are usually not equal to 
each other and their calculations are performed using Eq. (9), with the 
following relations of pressures and temperatures: 


for j 
if 


H3 : 


for j — HI 


PH3 1 

T b 

= ^H35 Pa = Ph 

Ph; 

T b 

= ^h; Pa — PH3, 

Ph; 

h 

= ^h; Pa = Phi, 

Phi; 

T b 

= Thi; p a = Ph 


( 22 ) 


The cross-sections, necessary for Eq. (9), are assumed as: 

^H 


7\hi = Au'i. 


m 




Pe = Pe^eR- 


(23) 
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and £hi = Ch 3 = 0.85. The signs of the mass streams rh H i and rh H 3 
will be defined later. 

The blower in the loop of hot air circulation assures that the 
Reynolds number calculated for the exchanger tubes is not less 
than 10 4 , when valves 1 and 3 remain closed. The Reynolds 
number caused by a flow from the first blower in such 
a situation is: 


Re H = 4m B1 /(7rd H MH n H) 5 (24) 

where rh B1 is calculated according to above-mentioned condition 
Re H = 10 4 . 

The Nusselt number for the proper heater tubes is calculated 
from the formula: 

Nu h = 0.022(Re H ) 0 ' 8 (Pr H ) 0 ' 6 - (25) 

where Ren is calculated as a varying value and uses rh H flowing 
through the heater. The value rh B1 is a small part of rh H , but exists 
when valves 1 and 3 are closed and the blower B1 works in the 
heater loop, Fig. 3. 

The coefficient of heat transfer equals to: 

^ah = Nu H A H /dH- (26) 

The heat stream delivered in the complete heater is a function of 
time and is determined as follows: 

0.AH = <*ah^h(?wh - Th)? (27) 


The essential equations for the proper part of the heater Vh are 
as follows. The equation of mass conservation: 


dp H _ 1 
dt V H 

where: 

• / if 

m H , { if 


and 



( ± ™H1 ± d7 H3 ), 


Phi > Ph, sign 
Ph > Phi , sign 


Ph>Ph 3 , sign 
Ph 3 >Ph, sign 



// // 


// // 



(32) 


(33) 


The equation of energy conservation has a form: 


cITh = 1 

dt Ph^h 


t„( ± m m ± iii„ 3 ) +«({} 


+m H 3^H3 I 0-AH 

-m H 3^H j J C v 


(34) 


The signs of the mass streams rh H i and rh H 3 are given in Eq. (33), 
but the term Qah is calculated for the complete heater. The pressure 
at the end of heater tubes is calculated as follows: 


where T w h - according to the described theoretical model - is the 
maximum temperature of the heater wall for a given rotational 
speed of the engine, and Th is a time-dependent temperature of the 
air inside heater tubes, at their outlet position. 

The pressures in both collectors are calculated in the following 
way. In the inlet, boosted collector H3, the pressure increases by 
ApH 3 , see - Eq. (31). The proper heater is composed of long tubes, 
which caused the pressure friction losses Apnf, calculated for the 
overall length of tubes. Therefore, the pressure in the outlet 
collector HI decreases by Apnf and is given in Eq. (39). 

Now, the essential equations for the heater will be presented. 

The calculation for the inlet collector Vh 3 is based on the 
essential equations for this part of the heater. The equation of mass 
conservation takes a form: 


Ph = Ph^h^- (35) 

Finally, the outlet collector of volume \4n is discussed. Its 

equation of mass conservation is: 

^- = y^(±m m (36) 

where: 

f if Ph > Phi , sign 
m Hl 1 „ „ 

l if Phi > Ph, sign 

and (37) 


dPH3 

dt 


—(±m 3 ±m H3 + m m ). 

V H3 


The signs of the mass streams are: 


m 3 


{ 


m H3 


{ 


if Pc>Ph 3, 
if Ph 3 > Pc, 

if Ph3 > Ph 
if Ph>Ph3, 


sign 

sign " —". 

sign // - // , 
sign " +". 


(28) 


(29) 


The equation of energy conservation is calculated from the 
formula: 

- ?H3 ( ± m 3 ± th H3 + m B i) 


dI H 3 _ 1 

dt Ph3^H3 


r if Phi > p E , sign // - // , 
l if Pe > Phi , sign " + 7/ . 


The equation of energy conservation takes a form: 


dT m = 1 

dt Phi^hi 


^hi (± /Phi =•= *n\ - tn B i) 


+ k(\ +rhm i' H } + ( + ™'1 E 

\l _m Hi^Hi J miTni 




(38) 


The signs of the mass stream rh H i and m j are given in Eq. (37). 
The pressure inside the collector volume I4u is calculated from 
the following equation: 


+ x 


+m 3 r c 1 

-rh 3 T H3 J 


+m H3 T H 

-rh H3 T H3 


+ di B iT H1 


(30) 


As the action of blowers is unidirectional, the positive mass flow 
rate of the blower B1 is assumed. 

Finally, the equation of the pressure pn 3 is: 


PH3 = PH3^H3^ + ^PH3- 


(31) 


Phi = PhiThiK - APHf- (39) 

At last, a cooler of the constant volume is presented. The cooler 
volume is divided into three parts, in a similar way to the heater: 

1. inlet collector -Vq\ 2 , 

2. proper cooler (a counter-current type) -Vci, 

3. outlet collector- V C 14 - 
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The volume of the cooler and the cooling surface are determined 
by: 

V a = ^cmci, (40) 

A a = 7rd cl / cl n cl . 

The working air is flowing from the expander E to the inlet 
collector Vcn, through valve 2 - see Fig. 4. Then it is transferred from 
the inlet collector through the proper cooler Vq to the outlet collector 
View, finally reaching the volume of the compressor C through valve 4. 

The blower B2 forces the working medium to transfer from the 
outlet collector to the inlet one in the cooler. We assume that the 
temperature in the outlet collector is transferred directly to the inlet 
and a pressure drop between them is covered by the B2 blower 
action. Additionally, the cooler collectors are considered to be settling 
chambers. The blower, which works in the cooling loop, causes that 
its Reynolds number is not less than 10 4 when valves 2 and 4 are 
closed. This number calculated for the cooler tubes being forced by 
the second blower (when valves 2 and 4 are closed) is determined as: 

Re cl = 4 rh B2 /(ird a p a ti a ), (41) 

where the value of rh B2 is calculated according to the above equa¬ 
tion for values of the Rea = 10 4 . 

The Nusselt number for the proper cooler part tubes is calcu¬ 
lated from: 

Nu a = 0.022(Re a ) 08 (Pr a ) 0 ' 6 , (42) 

where Rea is calculated as a varying values and uses rh a flowing 
through the cooler. The value rh B2 is a small part of rh a , but exists 
when valves 2 and 4 are closed and the blower B2 works in the 
cooler loop, Fig. 4. 

The coefficient of the heat transfer in the cooler is determined as: 

^aci = NiiciAci/da. (43) 

The heat stream in the complete cooler section is a following 
time-dependent function: 

0-AC1 = a AC\^C\(J w c\ - Tq), (44) 

where according to the described theoretical model, the tempera¬ 
ture T w a is the minimum temperature of the cooler wall, and Tq - 
of the air inside the cooler tubes, at their outlet. The temperature of 
the air inside the cooler falls, in fact, between values in the 
collectors C12 and C14, but the heat stream described by Eq. (44) is 
calculated for the overall length of its tubes. According to the 
theoretical model, the temperature is calculated at the outlet from 
the cooler tubes (the last plane of the proper part of the cooler) and 
is time-dependent only. 

The pressures in the collectors of the cooler are described 
similarly to the collectors of the heater. The inlet, boosted collector 
is now C12 and the outlet collector is C14. Similarly, the pressure 
increases in C12 by Apa2. but decreases in C14 by Apcif- 

The mass streams flowing in and out of the cooler are usually 
not equal to each other and their calculations are performed using 
Eq. (9), with the following relations: 


The approach applied to the cooler is identical to that one used 
for the heater and will not be repeated here. All the necessary 
parameters are determined from reasonable real values used in the 
simulation program. 


3.4. Engine power and efficiency 

The work of the expander and the compressor is calculated 
for the angle 0°<a < 360°, in the same way as for the 2-stroke 
engine: 

a = 360° 

L E = J p E (a)^da, 

a = 0° 

a = 360° 

L c = I p c (a)^da. 

a = 0° 

The above given pressures are determined while taking into 
account mass losses in the valves. 

The work of the engine is a sum of: 


(47) 


(48) 


L = L E + L c - L B i - ^B2 - U ) (49) 

where the values of work of the blowers are (i = 1,2): 

Tb i = 27rm B j (50) 


The symbols used are explained in Appendix 1, and an enthalpy 
isentropic increase of blowers is: 

h SBi = ^ 1 RT bBi (nj-l) (51) 


A complete power output of the externally heated 2-stroke 
engine is determined as: 



(52) 


The heat delivered inside the heater (following the theoretical 
model description) is: 


a = 360° 

Qh = / Qah^, (53) 

if the heat stream is calculated according to (27). 

The theoretical and practical efficiencies of the engine are 
equal to: 

V = L/Qa, 

(54) 

Vpr — ^/Ottpr- 


for j = C12 : 

if pa > Pci 2 : Pb = Pci 2 5 T b = Tq 2 ; Pa = Pcb 

if Pa > Pci 2 • Pb = Pch T b = ^C14; Pa = Pci 2 > ( 41:1 

for j = C/4 : ^ J 

if Pa > Pci 4 : Pb = Pci; T b = Tq] p a = Pci 4 ; 
if Pci4 > Pci : Pb = Pci4; T b = Tq 4 ; p a = Pa, 

The cross-sections and the coefficients necessary for Eq. (9) are: 
A:i 2 = A:i 4 = 7 r ^a n c//4; ?ci 2 = ?ci 4 = 0.85. (46) 


The efficiency may be defined in a different way. Calculation of 
the value of Qh is made without an additional heat exchanger. 
However, the temperature of the medium leaving the heater (Figs. 1 
and 6) is very high, so a pre-heater should be applied like that one 
described in [10]. The pre-heater rises the temperature of a gas in 
front of the inlet to the heater and causes considerable savings of 
the heating medium (decreasing the denominator in Eq. (54) from 
Qh to Onpr)- The efficiency of the EHVE, in this case, should increase 
significantly. 









2106 


J. Wojewoda, Z. Kazimierski / Energy 35 (2010) 2099-2108 


Table 1 

Rotation and thermodynamic parameters of EHVE. 


^ ehve [rpm] 

3000 

2500 

2000 

1500 

1000 

T wH max [K] 

1273 

1173 

1073 

973 

873 

Pm [kg/m 3 ] 

27.9 

28.0 

28.2 

28.3 

28.5 

Pc 14 [kg/m 3 ] 

14.0 

14.1 

14.3 

14.5 

14.9 


4. Results of the computer simulations 

The numerical task is reduced to a solution of the set of non¬ 
linear, ordinary differential equations (13), (15), (17), (19), (28), 
(30), (32), (34), (36), (38) with appropriate switching conditions, 
supplemented by a set of non-linear, algebraic equations. 
Stationary, periodic solution of this problem for engine operation at 
a given, constant angular velocity of the engine is presented. A 
solution of the custom written program built with the GCC 
compiler/linker is obtained using standard C++ language library 
procedures and solvers presented in [13]. The Runge-Kutta 4th 
order method is applied in the integration process. All simulations 
have started from identical, fixed initial conditions and continued 
until the solution reached integrated values differing in relative 
measures no more than 0.01% from those obtained at the end of the 
previous cycle. The numerical errors were within limits of accuracy 
of double precision numbers. Adaptive step-size integration 
method with forced decrease of step near switching values like 
valve opening/closing has been applied. The size of the integration 
step and limits of acceptance of the consecutive solutions has been 
chosen according to careful analysis of wide range of their values. It 
should be mentioned the system presents quick convergence of 
solution and rather small sensitivity to the step changes. Good 
agreement of the integrated values with experimental results 
published earlier in the paper [5] for the previous version of the 
engine makes the above statement satisfactory. Typically, stable 
solution has been reached after 100-150 cycles of simulated full 
crankshaft revolutions. The power of the described engine is 
calculated according to equations (49) and (52), the efficiency 
according to Eq. (54). The many of physical losses have been already 
taken into account. The heat transfer through cylinders and heat 
exchangers has not been considered. The losses caused by the heat 
transfer depend on a kind of insulation used and cannot be accu¬ 
rately defined at this stage. Therefore, it is estimated that total 
errors of the calculation in terms of effective power and efficiency 
can stay within 5% of results shown in Fig. 9. 

The calculations were performed for V<JV E = 0.75, considered as 
the optimal ratio for the investigated case when Vn ~ 5V E . The 
power and efficiency of the EHVE for this rate of the cylinder 
volumes are highest, but it is worth mentioning that both calcu¬ 
lated characteristics are flat in this region. Similar optimization and 
determination of Vq/V e is performed in [7], though physically both 
the blowers and their models introduced into the calculations are 
different. It is important that the above ratio is always less than 1 if 
we consider the engine with external heating. 

Other essential data were: 

d E - 0.083 m, s E = 0.0834 m, 

h 0E = 0.0025 m - dimensions of the cylinder E. 

d c = 0.083 m, s c = 0.06255 m, 

hoc = 0.0017 m - dimensions of the cylinder C. 

^lmax = ^3max = 0.001 m 2 , (tWO valves), 

^ 2 max = /Umax = 0.0005 m 2 - valves surface areas. 


Valves opening - index 0, closing - c, and respective transient 
angles - A: 


a 10 

- -10°, 

a ic 

o 

O 

00 

Aaq 

= 35 

^20 

= 160°, 

<*2c 

- 350°, 

Aa 2 

- 90 

^30 

- 110°, 

a 3c 

= 190°, 

Aa 3 

- 35 

&40 

- 210°, 

a Az 

- 360°, 

Aai4 

- 50 


The values of the above angles are established after preliminary 
optimization calculations. 

Dimensions of the heat exchangers (H - heater, Cl - cooler), 
d - pipe diameter, l - pipe length, n - number of pipes, V - volume: 

dj-j = 0.01 m, = 2.00 m, rij^ = 15, Vj-n = ^H 3 — 0.005 m^. 
d cl = 0.01 m, l C i = 2 . 00 m, n cl = 15, V cl2 = V C \ 4 = 0.005m 3 . 

A value of the Apn 3 depends on the design of collectors and heat 
exchangers. For simplicity, all values of the Ap’s are assumed to be 
10 kPa, when the EHVE works at 3000 rpm. 

77bz = 1.03 - input to output pressure ratio in the blower. 

TbBz - temperatures taken from the outlet of the collectors HI 
and C14. 

PBi = 0.7 - assumed efficiency value of the blowers. 

The temperature of the heater walls is assumed to decrease with 
the engine rotational speed and it is related to an amount of heating 
medium delivered, see Table 1. The minimum temperature of the 
cooler wall is constant and equal to T wC i = 300 K, based on standard 
parameters of cooling water supply. The coefficients of the heat 
conduction, dynamic viscosity and specific heats are taken from the 
tables of physical properties of the air. 

The electrical engines driving the blowers have to be controlled. 
Their temperatures and pressures at the inlet to the blowers, i.e., 
also the outlet collectors of the heat exchangers, decrease together 
with the rotation of the EHVE, but the density of the working air 
changes insensibly for different mass flow rates of blowers so the 
application of them is possible. The rotational speeds of the blowers 
are approximately constant and their angular velocities are 
cob i = 314 rad/s. The value of m B i > 0.06 kg/s is the minimum 
blower mass stream (of given heat exchangers) for which Re = 10 4 , 
when the valves of the engine are closed. 

The rotation of the EHVE, the maximum temperature of its 
heater walls, the averaged densities of working air in the collectors 
HI and C14 are shown in Table 1 for minimum mass flow rates of 
the blowers. 

The program based on the described model can be used for 
different calculations. Figs. 5-8 show some examples, explained in 
detail in their respective captions with the blowers for 
m Bi > 0.06 kg/s. 

The exception is Fig. 6 , which will be described here. According 
to the theoretical model we show here a temperature of the 
expander T E , a temperature of the compressor 7c, a temperature at 
the outlet from the heater Th, a temperature at the outlet from the 
cooler Tci, and temperatures in exchanger collectors - Thi, Th3, Tq2, 
Tq 4 . The accuracy of the applied approach to the heat exchangers is 
checked, e.g., the pressure and the temperature at the outlet from 
the heater H are always close the pressure and the temperature in 
the outlet collector HI. 

Fig. 9 shows power and efficiency of the engine with and 
without the blowers. The mass flow rate of both blowers are 0.1, 
0.2 kg/s. 

The heat delivered in presence of a pre-heater applied is smaller 
and is assumed to be QH P r = 0.85 Qh- This assumption is rather 
modest and can be easily met in a real installation. The efficiency of 
the EHVE in this case is obviously higher than in case without 
presence of a pre-heater and is presented in Fig. 9 as p pr . 
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5. Conclusions 


Appendix 1 


A two-stroke externally heated valve engine (called EHVE), 
entirely different from the Stirling engine, has been investigated 
numerically. Two counter-current heat exchangers, a heater and 
a cooler are applied in this engine. The engine may be supple¬ 
mented by two blowers (Figs, 1, 3 ,4), used to increase the 
volume of heat exchange when the valves of the EHVE are 
closed (Fig. 7). 

A lack of the blowers has been also considered. The EHVE has 
the moderate exchangers (up to 514), works in a closed Joule cycle 
with air as a working medium, and additionally has a conventional 
crankshaft and an oil lubrication system. 

A disadvantage of the EHVE realizing the Joule cycle is a short 
period of openings of valves Nos. 1 and 3, Fig. 7. It is a challenging 
but feasible engineering task. The statement is supported by the 
fact that intensive works connected with the valve timing are still 
developed, [15]. 

The calculations based on the model described are performed 
for V4/V e = 0.75, being the optimal value for which the EHVE power 
and efficiency are close to their optimum. Figs. 5-8 show pressures, 
temperatures, etc., described in their captions. 

The diagram given in Fig. 9 is essential for this paper. It shows 
the efficiency and the power of the EHVE if blowers exist or not in 
the installation as functions of rotational frequency of the engine. 

Taking into account that it is not an easy task to install the 
blower in the loop of high temperature and pressure, the profit of 
using blowers must be significant. In the present configuration, it is 
moderate, see Fig. 9. 

The efficiency of the EHVE does not increase practically. The 
power of the engine for m Bi = 0.1 increases by about 25%. Further 
growth in the value of m Bi has no practical sense, because the profit 
of using the blowers becomes insignificant. 

Nevertheless, the EHVE without the blowers reaches an 
efficiency level of the same value as the Stirling engine, when a pre¬ 
heater is applied, like in [10]. 

Additionally, the EHVE uses atmospheric air while the 
Stirling one works with helium as the working gas [10]. 
Considering this, the power produced by 1 L of the EHVE at 
1500 rpm, [10], when it works with the same gas, is calculated. 
The comparisons are following: the Stirling gives 51 kW, the 
EHVE without blowers - 59 kW. In this calculation, the 
conduction coefficients of helium are actually 5.3 times greater 
than these for the air. 

The engine may be applied close to garbage dams which 
produce bio-gas, for propulsion of small boats or yachts when they 
are supplied by tanks of bio-gases. It may be also used on small 
nuclear submarines. 

Figs. 5-8 show diagrams when the temperature of the heater 
varies in a range of 800-1000 °C. It is not always physically possible. 
For a solar heat source, so high temperature needs an application of 
a great set of special solar concentrators. The calculations of the 
engine for a lower temperature of the heater and rotational speed 
have been performed, but not published here. The engine operates 
properly in that case as well. 

Ships of the future will be propelled by artificial sails but this 
fails when there is no wind and the sun shines intensively. The 
described engine may be a unique source of power in such 
a case. 
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Notation 


A 
C v 
D, d 


h 

hj 

fioE* hoc 


hsBi 

L 

l 

mj 

P 

P 

Pr 

n 

d 

Qa 

R 

s 

Re 

t 

T 

V 

a 

a A 

a j,k 

\aj 

C 

V 

K 

X 

p 

OJ 


surface area [m 2 ] 

specific heat at the constant volume [J/kg I<] 

diameters [m] 

position of the piston [m] 

valve lift [m] 

clearances over the piston boundary positions for 

expander and compressor [m] 

entalphy isentropic increase of blowers [J/kg] 

mechanical work [J] 

length of tubes [m] 

mass flow rates [kg/s] 

pressure [N/m 2 , bar] 

power [W, kW] 

Prandtl number [-] 
number of tubes [-] 
heat [J] 

heat flow in the surface area [W] 
universal gas constant [J/kg K] 
piston stroke [m] 

Reynolds number [-] 
time [s] 

temperature [I<] 
volume [m 3 ] 

crankshaft angle [rad, grad] 
heat transfer coefficient [W/m 2 K] 
angle of valves [rad] 
increase in a valve angle [rad] 
loss coefficient [-] 
efficiency [-] 

isentropic expansion index [-] 
heat conductivity coefficient [W/mI<] 
density [kg/m 3 ] 

angular velocity of the crankshaft [rad/s] 


Subscripts 
a after 

Bl, B2 blowers 
b before 

C compressor 

Cl cooler 

C12, C14 collectors of the cooler 
E expander 

f friction 

H heater 

HI, H3 collector of the heater 
i number of blowers (=1,2) 

j number of a valve (=1, 2, 3, 4) 

k opening (/< = 0), closing (/< = c) 

max maximum 
w wall 
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